Induced long range dipole field enhanced antihydrogen formation 

in the p + Ps{n = 2) ^ e" + H{n < 2) reaction 



Chi Yu Hu, David Caballero and Zoltan Papp 
f— ^ ' Department of Physics and Astronomy, California State University, Long Beach, California 

O ■ 90840 

^ '. (February 2, 2008) 

Abstract 

CIh' We assume all interaction to be Coulombic and solve the modified Faddeev 

equation for energies between the Ps{n = 2) and H{n = 3), which involve six 

and eight open channels. We find that 99% of the antihydrogen are formed in 

j^ , H{n = 2). Just above the Ps{n = 2) threshold the S, P, and D partial waves 

^ [ contribute more than 4000 square Bohr radii near the maximum. Evidences 

^^ I indicate that the induced long range dipole potential from the degenerate 

^' Ps{n = 2) targets is responsible for such a large antihydrogen formation cross 

^ ■ section. 
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Due to the degeneracy of the excited hydrogen targets a long range dipole potential is 
induced in the field of the incoming charged particle. This potential is well known to be 
responsible for the formation of Feshbach resonances just below the n> 2 thresholds |1|J^. In 
the absence of relativistic effects and low lying open channels this attractive dipole potential 
can support an infinite set of bound states just below the threshold. Relativistic corrections 
remove the degeneracy and cut down the number of such states to only a few. The presence 
of low lying open channels embedding them in the continuum, they become resonant states 
and cause closely spaced oscillations in the cross section just below the threshold. 

Gilitis and Damburg [Q pointed out that due to the Levinson theorem ||^ similar oscil- 
lations exist just above the threshold. Our S-state cross sections in the energy gap between 
the Ps{n = 2) and H{n = 3) thresholds indeed reveal such oscillations and the antihydro- 
gen formation cross section is greatly enhanced in this region slightly above the Ps{n = 2) 
threshold. It is interesting to point out that the locations of the two Feshbach resonances 
below the H{n = 3) threshold coincide almost exactly to the two minimums of the anti- 
hydrogen formation cross section curve. In addition, antihydrogen formation cross sections 
for P and D partial waves are calculated for a number of energies. It is shown that the 
total antihydrogen formation cross sections from S, P and D partial waves rise almost up to 
14007raQ near the maximum region. 

The calculations are carried out using the modified Faddeev equations. The details of 
this method are given in Ref. [^]. We provide an outline below. The mass-scaled Jacobi 
vectors are defined as: 

Xa = r„(r^-r^), y„ = //„ r„ ^-^— ^^^ , (1) 

\ nifs + m^ J 

where Tq, = J2mf3m^/(mj3 + m^) and /ia = J2ma{l — nia/M), M = iria + rufs + m^, and 
(a,/3,7) are cychc permutations of (1,2,3). The mass of the antiproton used is mi = 
1836.1527, the mass of the electron and the positron are equal, m2 = m^ = 1 (the values are 
given in atomic units) . The Jacobi vectors of different channels are related by the orthogonal 
transformation: 
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^"~ [(M-m^)(M-m,)J ' ^^ 

Spa = {-lf~^sgn{a -(3){1- Cj^f/^ (4) 

In our particular case V3, the interaction between the p and e~, is a repulsive Coulomb 
potential which does not support two-body bound states. There are no asymptotic channels 
associated with this fragmentation. Consequently the total three-body wave function can 
be expressed in two components 

^ = V^i(xi,yi) +V^2(x2,y2). (5) 



Asymptotically ipi consists of an electron bound to the positron and a free antiproton, while 
i/j2, the rearrangement Faddeev component describes asymptotically an positron bound to 
the antiproton and a free electron. 

We use the modified Faddeev equations: 

-V^^%i^p,yp), (6) 

V^a = V3 + Vg , a 7^ /? = 1, 2, and Vi and V2, the interactions between the (e~,e+) and 
(p, e+) pairs, respectively, are separated into short- and long-range terms 

Vi''> = V^MUx^,y^), (7) 

V^'^ = V^iXa){l-Caix^,ya))- (8) 

The function Ca{xa,ya) vanishes asymptotically within the three-body sector, where Xa ~ 
ya ^ oo, and approaches one in the two-body cluster region, where x^ << ?/a ^ cxo. We 
use the following function having the required property: 

C(x, y) = 2{l + exp [{x/xoY/il + y/yoW' , (9) 

where u must be larger than 2 |^. In principle, Xq and yo are arbitrary, but they should be 
chosen to be consistent with the size of the scattering system for a rapid convergence. 
In bipolar basis ^ the Faddeev component is given in the form 

^a(x.,y„) = f: E E ^"'-^^-^"^ l^f^(x.,y.) (10) 

L=OM=-Lr+x=L ^^" 

where Yi{'^{xa,ya) = [y"'(^a) x ^r^'iVa)], a = 1,2, and 'ipa^^{xa,ya) is the partial compo- 
nent of the three-body wave function having total angular momentum L and relative angular 
momenta / and A associated with coordinates Xq and ja, respectively. Thus, for each L, the 
MFE 1^ is further reduced to a set of two-dimensional partial differential equations: 



-V^'\x^,y^) E {Yii''{£a,ya)\^^i^^^^^^^YA¥{x,,y,)), (11) 



where 
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Hl'^-Sl-Oiy-^^^, (12) 



and /3 7^ a = 1, 2. 



The summation with respect to angular momentum channels is truncated in Eq. (11). 
When the proper cut-off parameters Xq and t/Q in (9) is chosen the number of terms used for 
the S, P, and D partial waves ranges from 10 to 12 for the first Faddeev channel (p+Ps), 8 to 
10 for the second Faddeev channel (e~ + H). The optimal choices of parameters satisfy the 
conditions Xc,^^^ ^ IOxq; Va^a. ~ l%o, a = 1,2, where a;„„^^, ?/„„^^ are the respective cut- 
off distances which divide the asymptotic and interior regions. In the energy gap between 
Ps{n = 2) and H{n = 3) we used llOao and 75ao for Xq,„^^, a = 1, 2, and 450ao and 150ao 
foi'l/a^a.. respectively. 

In the MFE approach the total three-body wave function tjj^'^^ (x, y) , for each open channel 
0", is a vector having all ip^ as its components. We further split the wave function into 
interior and asymptotic parts using the vector equation 

ij^-\x,y)=F^-\x,y) + fi^'\x,y)+ ^ K^^Jj^sH^^v), (14) 

o-' = l 

where F^'^\x,y)\y^yy^^^^ = 0, f^^^\x^y) and fa'^'\x,y) are known incoming and outgoing 
asymptotic wave functions for the open channels a when y^ > ya^ax ^-nd they are equipped 
with spline continuity of values, first and second derivatives across the boundaries. They 
are otherwise identically zero in the interior regions. Kcr'a is related to the K-matrix by a 
kinematic factor 0. 

The partial cross sections between the incoming channel i and the outgoing channel j 
with total angular momentum L is given by 



a., = ^(2L+l) 



2K 
1 - iK, y 



(15) 



where ki is the momentum of the incoming channel. 

Upon substitution of (14) into (11), we obtain the vector equation 



^ma 



{H - E)F^-\x,y) = /. + E K^'J^'''\ (16) 

a'=l 

H is composed of two diagonal blocks of operators corresponding to the two Faddeev compo- 
nents on the left-hand side of (11) and two off-diagonal rectangular blocks involving the short- 
range potentials on the right-hand side of (11). la, l'^" ^ are known inhomogeneous column 
vectors. Replace the unknown vector F^'^\x,y) by F^'^\x,y) = U„{x^y) -\- YZ'P^l K„i^„U'^" ^ 
in (16). We find the f/'s satisfy the inhomogeneous equations 

{H - E)Ua = la (17) 

(ij_^)f/M=/W, a=l,...,a^,,. (18) 

The ?7's are solved using fifth order Hermite polynomial spline expansion 0-0. Matching the 
interior wave function F'^"\x^ y) with the known asymptotic wave functions at or near ya^r^ax 
we obtain highly over determined linear equations for the unknown Ka'a- Two independent 
methods are used for the solutions: (1) a standard least-square procedure and (2) projection 
on open channels. The details are given in [§,§|. 



The size of matrix equation used to solve (17,18) varies depending on the energy and 
the total angular momentum. They range from about 100000 x 100000 to 150000 x 150000 
for results accurate at least 5%. The solutions are obtained using a solver in ScaLAPACK 
library on the massive parallel Blue Horizon IBM computer at the San Diego Supercomputer 
Center. 

In the energy gap between Ps{n = 2) and H{n = 3) there are six open channels for S 
partial wave. The elastic cross sections are plotted in Fig. 0, the horizontal axis plots energies 
in Ps{n = 2) + p channel. The elastic cross section from e~ + H{ls) and p + Ps{ls) are 
relatively very small, they are not visible, that from e~ + H{2s) e~ + S{2p) are somewhat 
larger. They are plotted in Fig. |I], for comparison with the huge elastic cross sections 
from the channels p + Ps(2s) and p + Ps{2p). Due to the large induced dipole moment of 
the excited positronium, the Gilitis-Damburg oscillations are clearly seen in Fig. |1| below 
O.OOSRy, the energies are measured from the Ps{n = 2) threshold. The cross sections are 
rather monotonous between 0.003 — O.OOSRy. The structure below the H{n = 3) threshold is 
due to the two S-state resonances located at 0.8835Ry 0.8875Ry measured from the S{ls) 
threshold, or O.OOQORy and 0.0130Ry relative to Ps{n = 2) threshold. One notices that the 
oscillations of the two degenerate 2s and 2p channels are nearly 180° out of phase in the 
region of Feshbach resonances, although the effect is distorted by the statistical weight of 
1/3 on the 2p cross sections. 

Fig. D presents the S-state antihydrogen formation cross sections in the same energy gap. 
The formation cross section are greatly enhanced in the oscillation region below 0.003Ry. 
The largest cross section calculated is 2197rag. The two minima below the H{n = 3) threshold 
are located almost exactly at the location of the two Feshbach resonances ||,^ . 

Calculations of the P and D partial wave cross sections have been carried out for a 
number of energies in the oscillating region. There are 8 open channels and 64 partial cross 
sections, most of them will not be reported here. Table I presents all antihydrogen formation 
cross sections of S, P and D partial waves at four energies. There are 9 antihydrogen 
formation partial cross sections for S-wave and 16 such partial cross sections for P and D 
waves, respectively. We group them into three partial sums, a^^ , a}^ and cr]^ , where L is the 
angular momentum (S, P, D). o"^ is the sum of all partial cross sections with antihydrogen 
formed in Is state, a]^ is the sum of all partial cross sections having antihydrogen formed in 
2s and 2p states but originated from Ps(ls) target, a)^ is the sum of all partial cross sections 
with antihydrogen formed in 2s and 2p states but originated from excited Ps targets. 

Our calculations highlight three important points: 

1. The P and D partial waves follow the same tendency of strong oscillations near Ps{n = 
2) threshold similarly to that of S'-partial wave. Table I shows as as increases towards 
the Ps{n = 2) threshold, so do ap and a^. Seaton [|I| pointed out that the strength of 
the dipole potential from n = 2 target can support Feshbach resonances for all partial 
waves with L < 2. 

2. Since the cross section depend on the kinematic factor 1/kf in (15), where ki is the mo- 
mentum of the incoming channel, the oscillations are amplified near the lower threshold 
Ps{n = 2), and hence the antihydrogen formation cross sections are also greatly en- 
hanced. On the other hand, the Feshbach resonances just below the upper H[n = 3) 



threshold have neghgible contribution to the antihydrogen formation cross sections. 
In fact, Fig. ^ shows two minima in antihydrogen formation cross section located at 
0.8835Ry and 0.8875Ry above S{ls) threshold, respectively. Archer and Parker ^ 
found one 5'-state resonance at — 0.11492Ry below breakup, or 0.88454Ry above H{ls). 
Ho and Green ||^ found two S-state resonances at — 0.11606Ry and — 0.11206Ry below 
breakup, or 0.88340Ry and 0.88740Ry above H{ls). We used the threshold energy 
— 0.99946Ry for H{ls), due to finite p mass, to calculate all resonance positions. In 
comparison, our novel method to determine Feshbach resonance position seems to have 
at least three to four digits of accuracy! 

3. Not surprisingly, the long range polarization potential, mainly the dipole potential, 
couples strongly only amongst the excited 2s and 2p states which lie just below the 
targets Ps{n = 2). The ground states Ps(ls) and H{ls) contribute less than 2% of the 
antihydrogen formation cross section, and even less of them are formed in if (Is). As 
a result, near 99% of the antihydrogen are formed in 2s and 2p states with Ps{n = 2) 
targets. 

Archer and Parker |^ solved the Schrodinger equation using a form of hyperspherical 
coordinates. They calculated 5'-state partial cross sections up to H{n = 4) threshold. 
However, they purposely avoided the long range polarization potentials by projecting out 
only those with ground states targets and using cut-off hyperradius of p = 120. Oao. We 
used cut-off coordinates that are equivalent to p ~ 470ao. Our corresponding partial cross 
sections in the energy-gap between Ps{n = 2) and H{n = 3) are generally 25% lower. The 
complete coupling of all six open channels in our calculation may be responsible for the 
differences. 



Igarashi et. al. [^ made hyperspherical coupled-channel calculations for antihydrogen 
formation in a broad range of energies in p -|- Ps scattering including partial waves up to 
angular momentum L = 12 and included cross section from Born approximations for L > 13. 
Unfortunately there is only one of their energy points fell within the gap between Ps{n = 2) 
and H{n = 3) thresholds. This point is located at ~ 0.002Ry, its total formation cross 
section is approximately GAOna^, where our value at this point is approximately 3967raQ from 
three partial waves. This means that, at least at this energy, the lowest three partial waves 
contribute about 62% of the total formation cross section. Fig. and Table I show that this 
energy point is not in the region of optimal antihydrogen formation. Our largest calculated 
antihydrogen formation cross section is 13977rao from the lowest three partial waves. The 
energy at this point is 0.000494Ry. Table I shows that at this energy an /as = 3.15. This is 
less then the factor 2L + 1 in (15). But it is expected that higher partial waves contribute 
much less effectively at such low energy. The long range dipole potential from Ps{n = 2) 
targets is expected to cause oscillations for all partial waves with L < 2. We expect a much 
more drastic drop in (Tl/cjs for L > 3. Works are in progress to calculate contributions from 
L > 3. 
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TABLES 



TABLE L Antihydrogen formation cross 
p+Ps{n < 2) ^ e-+H{ls), p + Ps{n = 1) ^ 
The corresponding partial cross sections are 
L = S,P,D. The energy values (in Ry) are 
sections are given in vrag units. aL 
wave L. 



T.A' 



sections and partial cross sections for processes 
e~+H{n = 2) andp + Ps(n = 2) -^e-+H{n = 2). 
denoted by cf^ , aj^ and a\^ , respectively, and 
measured from the H{ls) threshold and the cross 
is the total formation cross section for the partial 
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FIG. 1. Elastic cross sections of the processes e +H(2s)^e +H{2p), p+Ps{2s) and p+Ps{2p). 
All energies are measured from the Ps(n = 2) threshold. The H{n = 3) threshold is located at 
0M395Ry. 
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Antihydrogen formation cross sections 

■ e~ + H{n < 2). The energies are measured from the 
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